Localized magnetic states in Rashba dots 
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We study the formation of local moments in quantum dots arising in quasi-one dimensional elec- 
tron wires due to localized spin-orbit (Rashba) interaction. Using an Anderson-like model to describe 
the occurrence of the magnetic moments in these Rashba dots, we calculate the local magnetization 
within the mean-field approximation. We find that the magnetization becomes a nontrivial func- 
tion of the Rashba coupling strength. We discuss both the equilibrium and nonequilibrium cases. 
Interestingly, we obtain a magnetic phase which is stable at large bias due to the Rashba interaction. 

PACS numbers: 73.23.-b,75.20.Hr,71.70.Ej 



I. INTRODUCTION 

Spin-related phenomena have recently attracted much 
attention, as they are the key ingredient in the new field 
known as spintronics 1 -. Two-dimensional (2D) semicon- 
ductors are appropriate materials to be used in spintron- 
ics applications since they offer the possibility of an elec- 
tric control of spins via tunable spin-orbit (SO) interac- 
tion. An important contribution to SO effects in 2D elec- 
tronic states of narrow gap semiconductors (e.g, In As) is 
the Rashba interaction 2 . This interaction is a generaliza- 
tion of the vacuum SO interaction from the Pauli equa- 
tion, H so — (eh 2 /4m 2 c 2 ) a ■ (VV(r) x k), which is small 
for non-relativistic momenta Kk -C mc, V(r) being the 
scalar potential. In semiconductors, the energy gap E g 
and the band splitting A are comparable in magnitude 
(Eg ~ A ~ leV) and, as a consequence, the SO coupling 
is enhanced by a factor m 2 c 2 /E g . 

The Rashba interaction is a type of SO interaction aris- 
ing when a 2D electron gas forms at the interface of a het- 
erostructure. To lowest order in momentum, the Rashba 
Hamiltonian reads 



H R = —([a,p y ] + a a 



[a,p x ] + a y ) , 



(1) 



where a is the Rashba coupling proportional to the elec- 
tric field producing the confinement. We take the con- 
finement direction along z. In Eq. {!])■ p = (p x ,p y ) is the 
2D momentum and <7, (i = x, y, z) are the Pauli matrices. 

We note that available experimental data 3 on few- 
electron quantum dots have been discussed in terms of 
Rashba spin-orbit coupling and exchange interaction 4 . 
Using the Spin Density Functional Theory it was showed^ 
that the competition of this coupling and the exchange 
interaction gives rise to the suppression of the Hund rule, 
and a dot with a closed configuration presents a param- 
agnetic behavior. We have to mention that these result 
have been obtained in the absence of the Coulomb inter- 
action. 

When the Rashba interaction is localized around a fi- 
nite region of a quasi one-dimensional ballistic wire [see 
Fig. [IJa)] i R ei - Ifl predicts the formation of quasi-bound 
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FIG. 1: (Color online), (a) The system under study consists 
of a quantum wire with a region of localized spin-orbit interac- 
tion of the Rashba type (a is the Rashba coupling strength) . 
Interactions are restricted to the Rashba dot. (b) Mapping of 
the upper system onto a quasi-localized level and a nonres- 
onant background channel with spin dependent couplings to 
external leads. 



states which are coupled to the nonresonant background 
channel. Both the potential well and the intersubband 
coupling is produced by the Rashba interaction alone. 
Furthermore, the quasi-bound states lead to enhanced 
backscattering, causing strong dips in the conductance 
curves of the wire as a function of the Fermi energy^ 
Since both the level position and broadening can be tuned 
with the Rashba strength a, these states are termed 
Rashba dots^. 

Recently, Lopez et al& have formulated a micro- 
scopic theory for transport across Rashba dots including 
Coulomb interactions in the dot. An important aspect of 
this model is that different regimes can be achieved by 
tuning the parameters of the Rashba Hamiltonian and 
this can be done modulating external electric fields ap- 
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plied to nearby gates. The difference between the Ander- 
son Hamiltonian 9 and the Hamiltonian proposed in Ref. 
H is twofold. First, in the Anderson Hamiltonian the spin 
is conserved at low temperatures, leading to the Kondo 
effect, but in Ref. [S| the Rashba dot Hamiltonian con- 
tains a spin-flip interaction because the localized states 
couple to the continuum states with opposite spins. Sec- 
ond, due to the Rashba induced precession term, the di- 
rect transmission channel presents a phase term similar 
to the Aharonov-Bohm case, but the phase is now spin- 
dependent^. Remarkably, despite these differences the 
system shows a persisting Kondo effect at low tempera- 
tures but with a novel gate dependence^. 

In this paper, we address the magnetic properties of 
Rashba dots. We follow Anderson's model for magnetic 
impurities in a metallic host and determine whether it 
is energetically favorable for the dot to form a localized 
magnetic moment. We show below that the Coulomb 
interaction can develop magnetic moments in a Rashba 
dot for a critical value of the ratio (Z7/r) cr it, which de- 
pends of the parameters of the Rashba interaction. Our 
results might also be important for quantum dots doped 
with magnetic impuritie o 1:L i 12 . Magnetic ordering in dots 
can be induced by the Coulomb interaction and the mag- 
netization can be electrically controlled even for a fixed 
number of electron a 13 ' 14 . 

This paper is organized as follows. We present in sec- 
tion II the model and calculate the Green function using 
the equation of motion method. The magnetic moment 
is determined in section III both for the equilibrium and 
the nonequilibrium cases. The main results are compared 
with exact numerical calculations in Sec. IV. The results 
are discussed in section V, which also contains our con- 
clusions. 



II. THEORETICAL MODEL 



We start with the model Hamiltonian: 
H = H e + Hd + H \y + H v 



vector k and spin a in the lead a = L, R. The nonreso- 
nant channel is described with the term Hw where the 
propagation phase acquired by a transmitted electron is 
spin-dependent (s CT = ±1 if a =T, !)■ Finally, localized 
and extended electronic states are coupled via the inter- 
action Hy. A pictorial representation of H is shown in 
Fig.ETb). 

The parameters of this Hamiltonian are U = U(a,l), 
V = V(a, I) and W — W(a,l), where I is the length 
of the Rashba induced square- well potential (we assume, 
for simplicity, that a(x) is constant if < x < I and 
zero otherwise)^, and cp = UrI with kn = ma /ft 2 . Im- 
portantly, these parameters can be externally controlled 
with gates by changing a and I. 

The form of H is similar to the Hamiltonian describ- 
ing the transport in a device formed by an Aharonov- 
Bohm interferometer with a quantum dot in one of its 
armsi^ but they differ in that the phases in the inter- 
action term W depend on the spin direction, and that 
each hopping process through the dot is associated with 
a spin- flip event. In the conventional Anderson model, 
spin is conserved, and this leads to Kondo correlations. 

In order to study the occurrence of magnetic mo- 
ments in this model we will calculate the Green function 
Gda{w) = ((do-lrfjj)), which obeys the equation, 



wG^H = l + (([d <T ,H]\dt)) 



(4) 



In the mean-field approximation, the spin-dependent en- 
ergy of the d-electrons is Hd,o = (d + U(n- a ). Using the 
general equation-of-motion method, we find 

£ *-f u - £fe 

ak 



AWcosip V 



S y-j-* u; - e k 



(5) 



(2) 



where S is given by the expression 



where 



OL,k,(7 



Hw = J2^ We ^ iPc U,a c R^ + H.c.) , 

k.a 

H V = J2(V<i )kitr d- a +H.c.). (3) 



a, k.a 



In this Hamiltonian we consider the spin quantization 
axis along the Rashba field (the y direction for transport 
along x), n a — d^da is the occupation number for elec- 
trons in the Rashba dot with spin a =f , J. and k is 
the creation operator of continuum electrons with wave 



s=1 - E 



w 



U - £fe 



(6) 



If we perform the summations over k in the relations 
above, the Green function becomes 

Gda{u) — [u - e da + «r + rV^cos^] , (7) 

where T = A/(l + x), x = -k 2 W 2 v 2 , A = irV 2 v, v being 
the continuum density of states at the Fermi level Ep (we 
take v as a constant function of energy). We note that 
the tunneling broadening A becomes renormalizcd into 
A/(l + x) due to the background channel when W ^ 0. 
Furthermore, the spin contribution due to the Rashba 
interaction is proportional to cos ip for the both spin ori- 
entations. 
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III. MAGNETIC MOMENT 
A. Equilibrium case 

The magnetization along the Rashba field direction is 
given by the difference between the occupancy expecta- 
tion value for spin up and spin down, 

m = (n dT ) - (n dl ) . (8) 

At zero temperature, the occupation reads, 



{n do 



p d (u))duj , 



(9) 



where pda{u) is the local density of states at the Rashba 
dot, defined in terms of the Green function as pda{oj) — 
-Im G d(T (w)/7r. 

Consider first the simple case U — 0. Then, 



^(r 2 + £M 2 )' 

where £ is given by 

= lu — €d + T^/x cos tp . 



(10) 



(11) 



Since U = 0, the energy £ is spin independent. Inserting 
Eq. (pi)]) in Eq. ©, we obtain 

/ \ 1 1 -i e d - E F - Y^fxcosLp 

(nd*) = - - - tan . (12) 

2 77 1 

Because (rid]-) = (ndi) even in the presence of spin- 
orbit interaction, we trivially have m = 0. As expected, 
equilibrium magnetic states arise due to the presence of 
Coulomb interactions only. 

Consider now the interacting case U ^ 0. We cal- 
culate the density of states using the spin dependence 
introduced by €da and we get: 



/a r 1 
PdA^) = -- 



(13) 



where 

C-o-(w) = w - Ed + Y^fx co$,(p - U{n-o) . (14) 
The occupation reads, 

1 1 _i e d - E F + U(n d -a) - Ty/x cos ip 



{n d a) = tan" 

2 77 



r 



(15) 

We analyze the formation of a magnetic state from the 
condition 



d{n dr7 ) 
d(n d -a) 



-Up d , a (E F ) 



(16) 



As a consequence, the condition for the magnetic state 
that d(n^} /d{n^ a ) < —1 becomes 



This relation is similar to the Stoner condition for the 
occurrence of the magnetic state in the itinerant-electron 
systems, and the correlations effects appear only as an en- 
ergy shift. A more accurate discussion, taking the energy 
dependence of A (via the density of states v) changes the 
magnetic region, which is known for a constant density of 
states. However, for our qualitative discussion we follow 
the wide-band approximation with a constant A. 

From Eqs. © and (fT5)) we find a pair of self-consistent 
equations for the magnetization m and the total electron 
density n d = (n dT ) + 



m = - > s (j cot , (18) 



rid 



^cot- 



i f (nd - s a m) - £(E F ) 



(19) 



From these two equations we calculate the size of the 
interaction above which a local moment develops. On 
the critical boundary describing the transition into the 
magnetic state, we approximately have m « and 
(nd]) ~ (ridi) ~ ri d - Thus, we find 



77 (1 



c 2 ) 



1 



(20) 



U P °d(E F ) > 1 



(17) 



where c = cot (7771^/2). This condition provides a num- 
ber of interesting predictions. First, for increasing x the 
function (C//A) cr i t decreases. Thus, the formation of 
magnetic moments is enhanced by the coupling to the 
continuum states, which is governed by the intensity of 
the Rashba interaction. Despite the fact that spin-orbit 
interactions are time-reversal symmetric and do not in- 
duce spontaneous magnetizations, indirectly the Rashba 
coupling makes it more favorable to generate magnetic 
solutions as compared to the case without Rashba inter- 
action. If x = (or, equivalently, W = 0) we recover the 
condition for the occurrence of the Anderson moments 9 . 
Second, (J7/A) cr ; t is a weakly function of the phase ip 
since in Eq. (I20|l the dependence on ip is only implicit 
through the total density rid- Nevertheless, in the gen- 
eral case the condition given by Eq. (I20p is far from being 
trivial since we recall that U , A and rid are complicated 
functions of the Rashba strength and the dot size&. 



B. Nonequilibrium case 

We now turn to the nonequilibrium case, where a finite 
dc bias V is applied between the two electrodes. The for- 
mation of magnetic moments within the Anderson Hamil- 
tonian out of equilibrium has been recently analyzed by 
Komnik and Gogolin, see Ref. [H. They find that the 
magnetic phase is stable at arbitrarily large voltages in 
the case of asymmetric couplings. Here, we assume sym- 
metric couplings (r^ = Fr) and show below that even 
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in this case the combination of Rashba interaction and 
finite bias leads to magnetic moment formation. 

At nonequilibrium, the spin-dependent occupations are 
given by the Keldysh (lesser) Green function G<(t, t') — 

i<4(tK(f)>, 



2^ 



(21) 



where G^(tv) is the Fourier transformed lesser Green 
function. We find 



r<l ^ _ o. ^sin^XJL - f R )/(l + x) + Cfc + f R )/2 

(22) 

where f a = 1/[1 + e™ u P"'] is the Fermi distribution 
function at the lead a = L, R with inverse temperature 
(3 = 1/kpT. Interestingly enough, the occupation de- 
pends on a term proportional to sin {utp) which has a 
different sign for opposite spins. This term appears only 
at nonequilibrium (fp ^ fp). Therefore, we expect a 
spin polarization (m 7^ 0) for a noninteracting Rashba 
dot (U = 0) induced by the interplay effect of external 
bias and Rashba interaction^ 

We take fi L = E F + eV/2 and /ip = Ep — eV/2 for 
the electrochemical potentials in the left and right con- 
tacts. As a result, we obtain a closed expression for the 
magnetization, 



1% 



■ sm tp 



tan 



.! eV - 2ti(E F ) 



2T 



- tan 



. x eV + 2£(E F ) 
2T 



(23) 



where T r = 4x/(l + x) 2 is the background channel trans- 
mission. We infer that the magnetization is negative for 
positive V, arising from the orientation of the effective 
Rashba field, which points along — y 18 ' 19 . The magneti- 
zation can be reversed if tp changes sign (equivalently, the 
Rashba intensity a). Obviously, the periodic dependence 
in Eq. (|23|) arises from the model but it is reasonable 
to assume a small a. Therefore, tp should not be very 
large. The periodic dependence on tp is obtained in anal- 
ogy with the Aharonov-Bohm effect and can be found in 
related spin-orbit systems (see, e.g., Ref. fl~7h . 
For tp = 7r/2 the minimum magnetization reads, 



= -2- 



JT r eV 
arctan — — 

7T 2T 



(24) 



which approaches —1 in the limits of T r — ► 1 and eV 3> 
2r. On the other hand, for gate voltages much larger 
than the applied bias the magnetization approaches zero 



1% TeV 



7T 



r 2 



(25) 



Finally, we note that the magnetization becomes finite 
and independent of the gate voltage in the limit of in- 
finite V , in which case we find the simple expression 
m = —\iT2. 



x=0, (|>=0 
x=0.1, (|>=0 
x=0.1,(|>=ji/4 _ 




FIG. 2: (Color online). Equilibrium magnetization of the 
Rashba dot for A = 0.2 and the values of x and tp shown in 
the figure. 



In the interacting case, one must replace £(w) with 
£- a (oj) in Eq. (|22[) . Then, the expression become in- 
volved and the full phase diagram can be obtained only 
numerically. However, the model is tractable in special 
cases. In particular, we focus on level energies around 
the particle-hole symmetric point (ed = —U/2). We de- 
fine the dimensionless parameters p — —Ed/U, y = U/T, 
z = eV/2T and R = 2{Y /U)\fx cos tp. Our goal is to char- 
acterize the critical line that separates the nonmagnetic 
and the magnetic phases. This is given by a curve p c 
versus y c in the p-y plane for different values of z. Then, 
for small values of p c — 1/2 we find (see the Appendix), 



(1 - Zz l 



7T 1 R 

2 {Pc -2 ) + TT^ 



(26) 

where Rq = \fx cos tp. This result shows that the phase 
diagram presents a dip for z > z* , where z* = l/v3 (as 
in Ref. [la ) but the form of the phase diagram is modi- 
fied by the Rashba parameter R. In fact, the dip position 
shifts away from the symmetric point due to the Rashba 
induced level renormalization. The most important con- 
sequence is that whereas m vanishes for Ed — —U/2 at 
large bias in the case without spin-orbit interactions^, 
the magnetization remains finite in the Rashba case. We 
have numerically checked this prediction (see below). 



IV. NUMERICAL RESULTS 

A. Equilibrium case 

We now numerically solve Eqs. {8} and (|15p . For sim- 
plicity, we neglect the dependence of the system parame- 
ters on a and treat U, A and Ed as independent constants. 
In Fig. [2] we show m as a function of Ed for Ep = and 
a fixed value of A. In the absence of spin-orbit inter- 
action (x = and tp = 0) the magnetization is zero for 
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FIG. 3: (Color online). Nonequilibrium magnetization of the 
Rashba dot for A = 0.2, x = 0.1, tp = tt/2 as a function of 
the level position for different values of the bias voltage V in 
the noninteracting case (U = 0). 



FIG. 4: (Color online). Nonequilibrium magnetization of the 
Rashba dot for A = 0.2, x — 0.1, ip = 7r/2 as a function of 
the level position for different values of the bias voltage V in 
the interacting case (U = 1). 



positive £d/U. When Ed/U decreases, there is a transi- 
tion point into the magnetic state, whose magnetization 
becomes maximal at the particle-hole symmetric point 
(sd = —17/2). At this point, m ss 0.64, in agreement 
with Ref. M- We now change the value of x and tp. These 
two parameters can be modified independently tuning a 
and I. Then, for nonzero x and tp = we find that 
the transition point shifts toward larger values of Ed/U . 
This results from the self-energy shift Ty/x cos ip found 
in Eq. ([7]). Moreover, we observe an increase in the am- 
plitude of to as x increases. This is a consequence of 
the Rashba coupling enhanced magnetic moment forma- 
tion discussed above [Eq. PD|) ]. Furthermore, keeping 
x constant and changing tp we find that the magnetiza- 
tion curve changes only slightly, confirming our earlier 
prediction. 



B. Nonequilibrium case 

The nonequilibrium magnetization for the noninteract- 
ing case (U = 0) is shown in Fig. [3] for increasing values 
of the external bias V. The curves are symmetric around 
Ed = 0, which corresponds to the alignment between Ed 
and the Fermi energy. The magnetization is nonzero for 
all finite values of V, as discussed after Eq. (|22|) . In the 
limit of large voltages, the curve becomes featureless ac- 
cording to Eq. (|2"3)). 

The interacting case is shown in Fig. [4j All energies 
are given in units of U = 1. For comparison, we also re- 
produce the curve corresponding to U = and eV — 0.1. 
At the same voltage in the interacting case, we observe 
that the magnetization curve follows the noninteracting 
curve for large Ed- This is reasonable since we are enter- 
ing the empty orbital regime for which interactions are 
unimportant. In the opposite regime, i.e., for negative Ed, 
strong correlations start to dominate and the interacting 
magnetization, although still finite, departs significantly 



from the noninteracting case. We obtain strong modi- 
fications in the magnetization curve for increasing volt- 
ages, favoring the development of magnetic moments due 
to the combined influence of interactions and spin-orbit 
coupling. For energies around the particle-hole symmet- 
ric point (Ed ~ —U/2), we find that the magnetization is 
reduced as V increases but, unlike the case without spin- 
orbit interactions, m does not vanish in the limit of large 
bias. This is excellent agreement with the prediction of 
Eq. (J26l). 



V. CONCLUSIONS 

We have studied the possibility of the occurrence of 
the magnetic models in a Rashba dot. The mean-field 
approximation has been used to calculate the magnetic 
moment and the critical value U/A for the occurrence 
of the local moments as function of the parameters of 
the Rashba Hamiltonian. This condition, expressed by 
Eq. (|20|) , is similar to the condition obtained for the An- 
derson model, 9 but contains also the parameter x which 
is determined by the Rashba interaction. Therefore, our 
calculation suggests a driving of magnetic moments by 
external electric fields via the Rashba interaction. We 
have demonstrated that the value of the local magneti- 
zation 777 at equilibrium depends on x, but it is worth 
noting that the curve m(ed/U) is not very sensitive to 
the change of the parameter <p. This result has been also 
shown in numerical simulations of the mean-field equa- 
tions. 

As in the standard Anderson calculation^ our mean- 
field approach breaks the local symmetry but in an exact 
solution accounting also for the effect of the spin fluc- 
tuations we should recover the spin rotational invari- 
ance. Nevertheless, even if the magnetic states found 
above are an artefact of the model, the mean field so- 
lution is interesting as such since it gives an indication 
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of the region of the coupling constants of the Hamilto- 
nian where the fluctuations give a relevant effect. Re- 
cently, magnetic moment formation was proposed as a 
mechanism to explain the temperature dependence of 
the conductance for different gate voltages in quantum 
point contacts , 20 ' 21 where the scaling behavior of the con- 
ductance close to pinch-off as a function of temperature 
was used as an argument for the Kondo effect occur- 
rence. Hence, our results can be useful for these systems 
when spin-orbit interactions become relevant. We believe 
that our calculation can also be important for magnetic 
semiconductors j 11 ' 22 . 

In the nonequilibrium regime, we have discussed the 
interplay between an external bias and the on-site inter- 
action energy when the spin-orbit interaction is present. 
The phase diagram we obtain is different from the 
nonequilibrium case studied in Ref. where the spin- 
orbit coupling was not considered. In Ref. HH, it is 
shown that the phase diagram presents a dip for z > . 
We have demonstrated that the spin-orbit interaction 
yields in Eq. (|26|) a correction given by the last term 
proportional to Ro- In the case of symmetric model 
(2e e i + U — 0) and large bias, the magnetization van- 



ishes in the absence of spin-orbit interaction*^ In con- 
trast, here we predict that the magnetization remains 
finite due to the Rashba interaction. Our numerical so- 
lution confirms this result, which can be particularly rel- 
evant for the experiments. It suggests that in materials 
with Rashba spin-orbit interaction the main contribution 
to the magnetization can be enhanced by applying a dc 
bias. 

As possible extensions of our model, an interesting pos- 
sibility is to take into account an energy dependent den- 
sity of states (specific for the semiconductors) like the 
gapless density of states p(u) ~ \lu\ t . This will give 
rise to an energy dependent A(w) and the resulting be- 
havior will likely differ from standard quantum dots^ 
Future investigations could also deal with the effect of 
correlations which was neglected in the present calcula- 
tions. Using the Hewson decoupling^ one might follow 
the method from Ref. [25| to calculate the effect of mag- 
netic correlations in the U — ► oo limit for systems with 
spin-orbit interaction. Finally, progress of experimental 
studies will be crucial for the directions development of 
this model. 



APPENDIX 



We present here a derivation of Eq. (f26|) . The mean occupations n j and have been calculated from the following 
relations: 



2im^ = — yT^sin^ 



2r 



2r 



7 r + tan-^±S(M-tan-- ey -^^ 



2r 



2F 



27m j = yTr-sincp 



taa _ 1 eV + 2 1 ,(E F ) +t&n _ 1 eV-^E F ) 



2T 



2T 



, eV + 2^(E F ) 
7T + tan -1 su ' - tan 

2r 



(A.l) 

. x eV-2^(E F ) 



2T 



(A.2) 



Using these results we obtain the magnetization m: 



1 r— 

m = — — \Tr sin if 
2n 



tan- zV + ^E,) tan _, eV-2 il{ E F ) ^ + ^ eV - 2^) 



2r 



2r 



2r 



2r 



i 

27 



tan- ^±MM _ tan- _ ^ eV + 2V^E F ) eV - 2^E F ) 



2T 



2T 



2T 



2F 



(A.3) 



From this expression we can see that for U — we have ^(E F ) = £,i(E F ) = £,{E F ) and the magnetization m has the 
value given in Eq. Eq. (T33]) with T r = 4x/(l + x) 2 . 

In the same way we calculate the total occupation n = + as 



n = — — v T r sin <p 
2n 



tan- £K±^M + tan -x *V-*l(E P ) _ tan _x eV + 2^E F ) _ eV - 2^E F ) 



2T 



2T 



2T 



21 



I 

~2~i 



tan- _ *Z^m + tan -: eV+^m _ ^ eV-2^E F ) ^ 



2T 



2T 



2T 



2T 



(A.4) 
(A.5) 
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which for {7 = becomes 



, 1 -l 8r ? 
n = 1 H tan — - r ; — 

7T 4r 2 - 4£ 2 + ( e y) 2 

for [(eV) 2 - 4£ 2 ]/4r 2 > -1. In the limit V -> this equation gives at £ F = and ip = tt/2, 



n = 1 tan 



T 2 -e 



2 ' 



(A.6) 



(A.7) 



From these relations we expect that the magnetic (n CT =/= n- a ) and non-magnetic solutions = n. CT ) exist for small 
V. In the following we will analyze the phase diagram taking into consideration the extra parameter V and the 
Rashba interaction. We introduce^^ the parameter n c , which runs from to 1, and the dimensionless parameters 
p = -e d /U,y = U/T,z = eV/2T,R = 2(T/U)y/xcosip. Thus, 



n c = 1 H — tan 1 [z + y c (p c + R/2 — n + n c )] tan 1 [z — y c (p c + R/2 — n + n c )\ 

7T 7T 



Derivating this expression with regard to n c we arrive at 

7T 1 



1 



y c 1 + [z + y c {p c + R/2 - n + n c )] 2 1 + [z - y c (p c + R/2 - n + n c )} 2 ' 

For p c = 1/2 we can fix n c = 1/2 and n = 1, 

7T _ 1 1 

^ ~ l + (z + i? ) 2 + l + (z-i? ) 2 ' 

where Rq = y^cosi^. We now write down the equation which contains the small parameter p c — 1/2, 

7T 1 1 



y c l + [y c ( Pc -l/2) + z + Ro} 2 1 + [y c (p c - 1/2) - z + i? ] 2 



Using Eq. (jA~T0| we can solve Eq. (jA.llj) iteratively, yieldin g Eq. (f26 



(A. 



(A.9) 



(A.10) 



(A.ll) 
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